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Abstract 

A new approach to the dynamics of relaxation and kinetics of thermahzation 
in a scalar field theory is presented that incorporates the relevant time scales 
through the resummation of hard thermal loops. An alternative derivation of 
the kinetic equations for the "quasiparticle" distribution functions is obtained 
that allows a clear understanding of the different "coarse graining" approx- 
imations usually involved in a kinetic description. This method leads to a 
systematic perturbative expansion to obtain the kinetic equations including 
hard-thermal loop resummation and to an improvement including renormal- 
ization, off-shell effects and contributions that change chemical equilibrium on 
short time scales. As a byproduct of these methods we establish the relation 
between the relaxation time scale in the linearized equation of motion of the 
quasiparticles and the thermahzation time scale of the quasiparticle distri- 
bution function in the "relaxation time approximation". Hard thermal loop 
resummation dramatically modifies the scattering rate for long wavelength 
modes as compared to the usual (semi) classical estimate. Relaxation and 
kinetics are studied both in the unbroken and broken symmetry phases of the 
theory. The broken symmetry phase also provides the setting to obtain the 
contribution to the kinetic equations from processes that involve decay of a 
heavy scalar into light scalar particles in the medium. 

ll.lO.-z; 05.30.-d; 24.60.-k 
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I. INTRODUCTION 



The motivation for understanding non-equilibrium phenomena in intermediate and high 
energy physics becomes more pressing with the experimental possibility of probing the 
quark-gluon and chiral phase transitions in upcoming heavy ion colliders at RHIC and 
large energy-luminosity hadron colliders at LHC In these high energy, high luminosity ex- 
periments (a/s > 200Gev/nucleon) a large energy density ~ IGev/fm^ is deposited in the 
collision region corresponding to temperatures T ?a 200Mev with time scales for the relevant 
relaxational processes of order t ^ Ifm/c. In these extreme situations, the challenge is to 
describe transport and relaxation phenomena from a fundamental microscopic theory since 
a description of non-equilibrium processes may lead to a better understanding of the experi- 
mental signatures of the transition and evolution of the plasma . Any model that attempts 
to describe the formation and evolution of the quark-gluon plasma and eventual hadroniza- 
tion in heavy ion collisions will have to succeed in providing a description of thermalization 
and relaxation during the expansion stage. The rate of thermalization and relaxation will 
determine if an equilibrium thermodynamical description is suitable or a full non-equilibrium 
treatment will ultimately be necessary to make quantitative statements on the evolution and 
signatures of the phase transition. Thermalization and relaxation are usually studied via 
a kinetic approach to non-equilibrium phenomena, typically through Boltzmann equations 
for the distribution functions. Such a description involves a wide separation of time and 
length scales and uncontrolled and drastic approximations are often made to simplify the 
problem. There is considerable effort in obtaining non-equilibrium equations that describe 
thermalization and relaxation starting from the microscopic theory that describes the full 
dynamics 0-0]. More recently this program has been extended to a detailed study of the 
kinetics of QCD and the excitations of non-abelian gauge theories [p HlO|] and the transport 
coefficients in scalar theories |jll] . 



Transport equations for quark matter are the subject of intense study since a deeper 
understanding of the thermalization and hadronic processes needs a full microscopic de- 
scription. A kinetic description based on QCD has been attempted over a decade ago but 
remains a formidable formal structure with very difficult implementation. While a descrip- 
tion from a fundamental theory is desirable, the technical obstacles are fairly formidable 
and instead a description in terms of low energy effective theories is emerging JT^|-[T^. 

A kinetic description of transport and relaxation is justified only when there is a wide 
separation between the microscopic time and length scales, namely the thermal (or Comp- 
ton) wavelength (mean separation of particles) and the relaxation scales (mean free path 
and relaxation time). The usual kinetic approach leading to Boltzmann transport equations 
involves Wigner transforms of two point (or higher point) correlation functions and even- 
tually a gradient expansion |]5| JT^jr3|| . Such a gradient expansion assumes that the center 
of mass Wigner variables are "slowly varying" but it is seldom clear at this level which are 
the fast and which are the slow scales involved. A "coarse graining" procedure is typically 
invoked that averages out microscopic scales in the kinetic description leading to irre- 
versible evolution in the resulting equations. Such averaging procedure is usually poorly 
understood and justified a posteriori. 

It is recognized that at temperatures much larger than the masses of the particles, it is 
necessary to perform a resummation of the perturbative series to incorporate the relevant 
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microscopic time scales in the description of thermal processes. Such a non-perturbative 
resummation scheme has been proposed by Braaten and Pisarski and since then used 
to obtain damping rates [0,|1^]. A program to incorporate self-consistently the microscopic 
time scales within a kinetic approach has been proposed by Lawrie in a manner that is 
similar in some respects to the hard-thermal loop resummation. In the quark-gluon plasma 
this program acquires further significance, since for temperatures T ^ 200Mev the thermal 
masses of the "light quarks" are of order Mq{T) ^ gT > rriq ^ 5 — lOMev with rriq the 
"current" quark masses and g the quark-gluon coupling constant P JlO|JT5 |. 

Another realm in which relaxation and thermalization also play a fundamental role is the 



description of the reheating stage in cosmological inflationary scenarios [p!q-pO|. Recently 
new non-equilibrium phenomena have been recognized to play a role during the reheating 
stage in the post-inflationary epoch in early universe cosmology. These phenomena are 
the result of profuse particle production during the stage of parametric amplification of 
quantum fluctuations for large amplitudes of the inflaton field PT|-Z1|. A similar mechanism 
has been recently proposed to play a role in hadronization during the supercooling stage of 
the quark-gluon phase transition [E^ . Since the particles produced during this stage are in a 



non-equilibrium distribution the final stage of reheating needs a deeper understanding 
of thermalization processes. 

The reliability of any kinetic approach to the description of thermalization and relaxation, 
either in early universe processes or hadronization and equilibration during the quark-gluon 
plasma transition, hinges upon a complete understanding of the time scales involved. Since 
in extreme environments the microscopic time scales are modified by the medium, it is 
important to incorporate these effects in the kinetic description. 

The purpose of this article is to provide an alternative program for the description of 
relaxation and kinetics in scalar field theories. The relevant microscopic and relaxational 
time scales are incorporated in a consistent first principles calculation starting from the 
microscopic theory. Such an approach allows us to identify the several different "coarse 
graining" procedures with a detailed understanding of the relevant scales. The resummation 
of the hard thermal loops is incorporated consistently in the derivation of the equations of 
motion that describe relaxation and kinetics and provides a systematic scheme to include 
off-shell and higher order corrections to the resulting kinetic equations, including terms 
that provide departure from chemical equilibrium on short time scales, and renormalization 
effects. We find that the hard thermal loop resummation modifies dramatically the estimate 
of the scattering rate in the Boltzmann equation for long wavelength modes, from that of a 
"naive" or (semi) classical argument. 

Furthermore, our approach provides a direct proof of the relation between the relaxation 
processes of quasi-particles and the thermalization rate of the quasiparticle distribution 
function in the "relaxation time approximation" of the Boltzmann equation that describes 
their evolution. Such a relation was proposed by Weldon for the identification of the quasi- 
particle damping rate ^6] but to our knowledge a complete proof of such relation has not 
been given within the context of the Boltzmann equation as derived from a microscopic 
theory. 

Although we focus here on a scalar field theory both in the broken and unbroken sym- 
metry phase to distinguish the different scales and processes, the methods can be extended 
to other models and theories almost straightforwardly (obviously with the complications of 
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gauge invariance in gauge theories). The broken symmetry case allows us to study the situ- 
ation in which the scalar field is coupled to other lighter scalars with a trilinear coupling and 
to include decay and recombination processes in both the relaxation and kinetic description. 

In the next section we introduce the model and the techniques. In section III we study the 
relaxation of "quasi-particles" obtaining the real-time equation of motion in the linearized 
approximation, including consistently the hard-thermal-loop corrections to the self-energy. 
In section IV we introduce a new formulation of kinetic theory that incorporates the hard- 
thermal-loop corrections consistently, displays all the relevant time scales and leads to a 
well defined "coarse-graining" procedure. This formulation allows a systematic improve- 
ment on the kinetic description and, if necessary, consistently includes off-shell effects and 
renormalization in a computationally accesible form. 

Section V presents our conclusions, summarizes the main results and poses new questions. 



II. THE MODEL AND TECHNIQUES 

We will restrict our study to a self-interacting scalar theory both in the broken and 
unbroken symmetry state. The Lagrangian density is given by: 

^0 = \ {d,^f - - . (2.1) 

with mo ; Aq the bare parameters. 

As mentioned in the introduction, the first step towards understanding the kinetic regime 
is the identification of the microscopic time scales in the problem. In a medium, the bare 
particles are dressed by the interactions becoming "quasiparticles" . One is interested in de- 
scribing the relaxation of these quasiparticles. Thus the important microscopic time scales 
are those associated with the quasiparticles and not the bare particles. If a Boltzmann 
equation is obtained in some perturbative scheme, such a scheme should be in terms of the 
quasiparticles, which already implies a resummation of the perturbative expansion. This is 
precisely the rationale behind the resummation of the hard thermal loops in finite tempera- 
ture field theory and also behind the self-consistent treatment proposed by Lawrie [jl7 



In this scalar theory such a resummation can be implemented by writing in the Lagrangian 

ml = m\{T) + 5m^{T) (2.2) 

where m\{T) is the renormalized and temperature c?ej»enc?ent quasiparticle mass which enters 
in the propagators and 5m?{T) is a counterterm which will cancel a subset of Feynman 
diagrams in the perturbative expansion and is considered part of the interaction Lagrangian. 
This method is the simplest form of implementing the hard thermal loop resummation in the 
scalar case and was used within this context by Parwani . We note that this renormalized. 



temperature dependent mass determines the important time scales in the medium but is not 
the position of the quasiparticle pole (or, stricly speaking, resonance: see the discussion 
below). We could also introduce counterterms for wave-function and coupling constant 
renormalization and proceed to a perturbative expansion of the BPHZ type in terms of 
the renormalized couplings and fields. Such a possibility will be discussed later within the 
context of improvements on the method and renormalization. 
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A. Non- Equilibrium Techniques 



The field theoretical methods to describe processes out of equilibrium are known and 



described at length in the literature [p^-|3^ . The basic ingredient is the time evolution of an 
initially prepared density matrix, which leads to the generating functional of non-equilibrium 
Green's functions in terms of a path integral representation along a contour in the complex 
time plane. This contour involves a forward time branch, a backward time branch and a 
third branch down the imaginary time axis to time t = —if3 if the initial density matrix 
describes an equilibrium ensemble at initial temperature 1/(3. For the computation of real- 
time correlation functions at the order of approximation we consider, the only role of the 
imaginary time branch is to determine the boundary conditions on the propagators. 

The fields living on the forward and backward branches will be labelled with + and — , 
respectively, and the effective Lagrangian that enters in the path integral representation of 
the non-equilibrium generating functional is given by 

£none^ = C[<t>^]-C[<^-]. (2.3) 

From this path integral representation it is possible to construct a perturbative expansion 
of the non-equilibrium Green's functions in terms of modified Feynman rules. Correlation 
functions are obtained as functional derivatives with respect to sources on the respective 
branches. The non-equilibrium Feynman rules are: 

i) The number of vertices is doubled: those in which all the fields are on the -|- branch 
are the usual interaction vertices, while those in which the fields are on the — branch have 
the opposite sign. 

ii) The combinatoric factors are the same as in usual field theory. 

iii) The spatial Fourier transform of the (bosonic) propagators are 

= G>it,t'Mt-t')+G<it,t')eit'-t) , (2.4) 

= G>{t,t'Mt'-t)+G<{t,t')Q{t-t') , (2.5) 

G^M = -G<{t,t') , (2.6) 

G,+{t,t') = -G>{t,t'), (2.7) 

G>{t, t')=ij d^x e"*^-^ (<l>(f, t)$(0, t')) , (2.8) 
G<(t,t') = 1 fd^x e-*^" {<l>{6,t')<l>{x,t)) . (2.9) 



where $ is the bose field. Now we have to specify the properties of the initial state. Par- 
ticularly convenient is the choice of a thermal initial state with temperature T. Then the 
density matrix of this initial state is p = e~^'"''^ , where Hq is the Hamiltonian for times 
t < 0. This choice of the initial state determines the boundary conditions on the Green's 
functions. These are the usual periodicity conditions in imaginary time (KMS conditions): 

G< (f , t; x', t') = G> (f, t - x', t') (2. 10) 

Finally, the free-field Green's functions for bose fields are constructed from the following 
ingredients: 
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GUt,t') = -^{[l + n(^,)]e-^-'=(*-*')+n(a;fc)e^-^(*-*')} , (2.11) 



2uJk 

GUt,t') = -^{[l + n(^,)]e^-'=(*-*')+nK)e-^-'=(*-*')} , (2.12) 



Uk = \Jk'^ + m\{T) , nb{u}k) = ^^^^ _ ^ , (2.13) 

An important property that will be used in the calculations that follow is the relation 

G>{t,t')=G<{t',t) (2.14) 

III. LINEAR RELAXATION 

In this section we study the relaxation of quasiparticles, from the linearized real-time 
equation of motion for the expectation value of the scalar field 

(P{x,t) = {<^+ix,t)) = {<l>-{x,t)) (3.1) 

where <l>^ are the fields on the forward and backward branches of the contour respectively 
and the expectation value is in the ensemble of the initial density matrix. We derive the 
equation of motion for 0(x, t) using the tadpole method [^,0], which consists in splitting 
the field $ into its expectation value and a fluctuation as 

$±(x,t) =0(x,t) +^±(x,t) , (3.2) 

where are the fluctuation operators defined along the respective branches. The effective 
equation of motion for the field (j){x, t) follows from the condition: 

<i^^{x,t) >=0 . (3.3) 



This tadpole method is different from that used by Elmfors et al |]3^ in that the present 
method is implemented with the non-equilibrium, real time functional and Green's functions. 
After the shift (|3.2|), the non-equilibrium Lagrangian density reads 



-^W(T)(20^+ + (^+)2)} - ^ (3.4) 

with /^ofV'^] the free field Lagrangian density of a field with mass m^(T). To implement 
the tadpole method one considers the linear, cubic and quartic terms and counterterm as 
perturbations and the condition ( |3.3| ) is imposed order by order in a perturbative expansion, 
using the Feynman rules given above p3|j3^ 
We now analyze different cases in turn. 
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A. T >> Tc'. Unbroken symmetry case 



We are interested in studying relaxation and kinetics both in the disordered high tem- 
perature phase as well as in the ordered low temperature phase. This is achieved by allowing 
the zero temperature renormalized squared mass to be negative, that is 

ml{T = 0) = -ImliT = 0)\. (3.5) 

When the temperature is much larger than the renormalized (zero temperature) mass, the 
hard-thermal loop resummation is needed to incorporate the physically relevant time and 
length scales in the perturbative expansion. In the equation of motion this is achieved by 
requiring that the mass counterterm 5m? {T) cancel the 1-loop tadpole contribution. 
This leads to the following self-consistent "gap" equation for m\{T) p7| , p6| : 



ml{T) = -ml + - [ ^ , |l + ^ ^ 1 (3.6) 



The divergence in the zero temperature part of (|3.6| ) (quadratic and logarithmic in terms 
of a spatial momentum cutoff) can be absorbed in a renormalization of the bare mass by a 
subtraction at some renormalization scale: 

2 ,m\ I 1 I f dk k^ ( 2 1 ^ , , 

mp T = — mp H — / — - — , = < 1 H , > — subtraction (3.7) 



where both and the subtraction depend on the renormalization scale and the renormal- 
ization prescription. Now the gap equation is finite. For T ^ m/j(T) we obtain 

^liT) = -\ml\ + ^ 



AT 



+ O (Air") for T > v^T > |m/j| (3.8) 



24 

The leading term of the last expression of (|3.8|) provides the correct microscopic time scale. 



In the massless case it serves as an infrared cutoff for the loop integrals p7| , p5| . We will 
concentrate on the linear relaxation case in which the equations of motion for the expecta- 
tion value are linearized in the amplitude of (f){x,t). This approximation is valid for small 
amplitude fluctuations from the minimum of the finite temperature effective potential in 
the unbroken symmetry phase. Two-loop diagrams contributing to the equation of motion 
obtained from the condition = are shown in figure 1. This condition leads to 

4){x\ t') - V^0(f' , t') + m^(T)0(f' , t') + / £x"dt"^ret{x' - x\ t' - t")(j){x", t") = (3.9) 



The mass counterterm has been used to cancel the 1-loop tadpoles. To two-loop order we 
find 

J:ret{x- x',t-t') = -i—\{-iG>{x-x',t-t')f-{-iG<{x-x',t-t')f]e{t-t') (3.10) 

6 
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The equation obtained from {ip~) = is the same as ( p.9|) as a consequence of unitarity. 

Translational invariance of the self-energy makes it convenient to write the equation of 
motion for the spatial Fourier transform of 0. We therefore introduce 

4>{x.t) = J Mt) (3.11) 

Eretix -x,t-t') = J ^ e^^-^-"^'^ Epit - t') (3.12) 
In order to solve the integro-differential equation, we will impose the initial condition 



(p^{t < 0) = such that this configuration is "released" at time t = [Q. Under this 
condition the evolution equation becomes 

0p^(t) + cu20p.(t) + f dt'll^it - t')Mt') = (3. 13) 



where 

= \P\' + ml{T) (3.14) 

T.p{t -t') = -- -— I i27rf5'ip- k,-k,- k,) 

{sin [{uki + uJk2 + ^ks){t - t')] [(1 + nkj (1 + rik^) (1 + Uk^) - Uk-, Uk^ Uk^ ] 
+ 3 sin [(cjfc, + ujk2 - ^ki){t - t')] [(1 + ^fci) (1 + ^fca) - n-k, Uk^ (1 + ^^3)]} 

(3.15) 



The linearized equation of motion ( p.l3| ) can now be solved by the Laplace transform. With 



the boundary conditions (j)^{t = 0) = (pi-^; (f)^{t = 0) = and denoting the Laplace transforms 
of 0p*(t) , ^p{t) by 0p-(s) , Sp>(s) respectively (here s is the Laplace transform variable) we find 

The time evolution is found by inverse Laplace transform, i.e. by integration in the complex 
s-plane along the Bromwich contour 

/ico+e /7q 
-ioo+e Zm 

with e to the right of the real part of all the singularities of the Laplace transform 0p(s). 
The Laplace transform of the retarded self-energy has the following spectral representation 

^p(^) = " / :^-r^2 dpo . (3.18) 

where p{po, p; T) is the spectral density, given to two- loop order by 
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p{k,po;T) = — / — — — — — — {2ny6%k-ki-k2-h 

6 J (ivr)'^ (iTr)'^ (iTr)'^ 8c<Jic<j2t^3 

{5(po - t^fei - t^fc2 - uJks) ((1 + r^fei) (1 + nk2) (1 + rifcs) - nki nk2nk^) 

+3 5(po + t^fei - t^fe2 - t^fcs) (^fci (1 + nk2) (1 + n^s) - (1 + nfci) ra^a ^fcs)} 

(3.19) 



This spectral density displays two different type of contributions P6| , p!6| , |37| : the first one 
(corresponding to the first 5-function) corresponds to the production process (fxjxf) with 
the Boltzmann weight (1 + nki){l + nfc2)(l + Uks) minus the recombination 000 ^ in the 
medium with the weight nkink2nk^. The second contribution (second ^-function) corresponds 
to the scattering process 00 00 with weight nki{l+nk2){)-+nkz) minus the inverse process 
with weight nkink2{^+nk^) . The self-energy Sp-(s) has an imaginary part along the imaginary 
axis in the s-plane which can be read off from the spectral representation: 

S/,p>(iu; ± 0+) = ±7r sign (t^) \p{\uj\,p,T) - p{-\uj\,p,T)] (3.20) 

where the sign(ci;) reflects the retarded nature [^,0. To two-loop order we find the imagi- 
nary part to be given by: 



y.[5{\uj\ - iOki - ujk2 - i^ksWl 

+36{\uj\ + UJkl - UJk2 - ^k-i)M2 

+5{-\uj\ - Ukl - UJk2 - UJksW^ 

+6{-\u\ + Uki - ujk2 - uJk-sWi] (3.21) 

where 

M = (1 + Uki) (1 + nk2) (1 + riks) - riki nk2 riks 
J^2 = riki (1 + rifca) (1 + ^fcs) - (1 + «fci) ^^^2 nks 
A/'a = Uki nk2 riks - (1 + ^fci) (1 + ^^2) (1 + riks) 

Xi = {1 + Uki) nk2 riks - TLki (1 + nk2) (1 + nkz) (3.22) 

This result coincides with that found by Wang and Heinz [^]. The first term gives a 
contribution above the three particle threshold at > 3mij(T), the third term vanishes 
identically, the second and fourth term give the coUisional contribution and have support 
all along the imaginary axis [H^. Therefore along the imaginary axis we find 



d^ki d?k2 d^k. 



S,^(^a; ± 0+) = S^^g(.;) ± ^S,^g(.;) (3.23) 

The real part can be obtained by a dispersion relation. It is very difficult to obtain in 
general, but we find that it is quadratically and logarithimically divergent. These divergences 
can be subtracted away from the dispersion integral and the subtraction absorbed in a mass 
and wave function renormalization. 

In order to invert the Laplace transform and obtain the real-time evolution, we need to 
understand the analytic structure of (f)^{s) in the complex s-plane. Isolated single particle 
poles would be at the values s = ±iQ{p,T) with Q{p,T) the solutions of 
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- n^p, T) + \p\^ + Ml{T) + S(s = ±1^1, T) = 



(3.24) 



Since the self-energy is complex along the imaginary axis, such a pole solution would be 
complex. However it is easy to see that there is no solution to ( p.24| ) in the first (phys- 
ical) Riemann sheet. The pole moves off into the second (unphysical) Riemann sheet as 
corresponds to a "resonance" (or quasiparticle pole) |^. The discontinuity of the Laplace 
transform across the imaginary axis is given by: 



0(P) s = iui + 0^) — 0(p, s 



lUJ 



0^ 



;.p- 2iujT.j{uj,p,T) 



\p\^ - Ml - T.R{uj,p,T)f + ^i{uj,p,Ty 



(3.25) 



This discontinuity is the only singularity in the (first Riemann sheet) complex s-plane and 
it vanishes at the origin. This fact allows us to deform the integration countour to wrap 
around the cut. 

Finally for the time evolution we obtain 



0p(t) 



a;S/(u;, p; T) cos(u;t) duj 



SK(a;,p;T)]2 + S,(a;,p;T)2 



(3.26) 



with c<jp given by ( |3.14| ). 

Setting t = in the above expression leads to the sum rule 



2 

TT Jo 



ijjY,i{ijJ^P\ T) duj 



(3.27) 



which has been previously obtained by Pisarski I^Sf in a rather different manner. 

For very weak coupling the integrand in ( |3.26|) features a narrow resonance at u"^ = 
ujp + (9(A^), and can be approximated by a Breit-Wigner form leading to the real-time 
evolution [|33| 



with 



0p-(t) = 0,;p-Z[p;T]e-r®^)* 
T,j{uj = ujp,p;T) 



cos{ujpt + a) 



T{p;T) ^ Z[p,T]- 



Z[p-T] 



dJ:n{uj,p,T) 



a 



dT.i{u,p,T) 



1 + 0{\^ 



(3.28) 

(3.29) 
(3.30) 
(3.31) 



where Z\p\T\ is the wave function renormalization defined on shell. To the order that we 
are working, this will be set to one in ( p.29| ). r(p; T) is the coUisional relaxation rate. The 
Breit-Wigner approximation, however, is valid only for times t < ^ In (^^). For longer times 
the fall-off is with a power law determined by the behavior of the spectral density at small 
frequency ||33[| . 
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The calculation of the decay rate T{p; T) is in general very difficult for arbitrary momen- 
tum and temperature (see however |^). For p -C Mr{T) -C T we can approximate it by 
the zero momentum limit, which can be calculated relatively easily. The only contribution 
to the imaginary part on-shell is given by the second and fourth terms (coUisional terms) in 

(Q). 

We find for T > ^AT > |m/j|: 



r(p=0;T) 



64V24 



(3.32) 



The same result had been obtained previously by many authors [27,35,37 



B. T « Tc'. Broken symmetry case 

In this section we consider the linear relaxation of fluctuations around the broken sym- 
metry state, for non-zero temperature but below the critical value. In this case the scalar 
field acquires an expectation value v which at tree level is given by 




(3.33) 



The small amplitude fluctuations around the broken symmetry state have a tree level mass 
squared given by 

Ml = ml-^ = 2\ml\ (3.34) 

We now write the field $^ as 

$±(f,t) = t, + 0(x,t) +V^±(f,t) . (3.35) 

We will not fix v at the tree level value given by eq. |3.33| , but instead the v.e.v will be 
determined order by order in the perturbative expansion by requiring that the linearized 
equation of motion for 0(x, t) be homogeneous. Furthermore as in the previous case we will 
write Ml = Ml{T)+5m^{T) and treat 5m^{T) clS cl counterterm. The microscopic time scale 
is determined by Mj^^(T). The condition {'tp'^{x,t)) = implemented order by order leads 
to the equation of motion for 4>{x,t) as before [^. In the broken symmetry phase there 
are new Feynman diagrams that contribute an absorptive part to the self-energy at one- 
loop order which are proportional to the v.e.v. Now there are new vertices in the "shifted" 
Lagrangian density. Those that contribute to the linearized equation of motion are obtained 
from the part of the Lagrangian density given by 

= !^vi;^ml - ^v') - y (30(^+)^ + (V^+)3)| - ^ } (3.36) 

To one-loop order there are two tadpole terms that contribute to the equation of mo- 
tion. One is independent of and is absorbed in a renormalization of v (since it gives an 
inhomogeneous contribution to the equation of motion). The other is proportional to and 
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is cancelled by the mass counterterm leading to the "gap equation" ( p.6| , p.7| ), but in the 
present case Mn{T) > T > Vat. The remaining one-loop diagrams contributing to the 
self-energy with an absorptive part are shown at the top of figure 2. We obtain the equation 
of motion given by eq. ( p.9|) with a new one-loop contribution to the retarded self-energy as 



well as the two-loop contribution obtained in the previous section and given by eq. (p. 101) 



We find the one loop contribution to be given by (see also [p3l , p9| ) 

AV 



-li-tG^ix -x',t- t')y - {-iG<{x -x',t- t')y]e{t - 1') 



(3.37) 



with the following spatial Fourier transform: 



2 J {27if2uj,uj^,_^^ 
+2nk sm[{ujj:^j.- ujk)it - t')]} 



{(1 + 2nk) sm[{u^^^+ujk){t - t')] 



(3.38) 



The Laplace transform can again be written in a spectral representation as given by 
( p.l8| ), and we find the one- loop contribution to the spectral density to be 



{2tiY2uj^uj^ 



{(1 + 2nk)5{po - ^fc - ^k+p) - '^^k^iVo -^k + ^k+v)} 



k'^k+p 



k+p' 



(3.39) 



The imaginary part is obtained from ( p.20|) , and after analyzing the kinematical regions we 
find for a; > (for a; < there is a sign change) (see also p9|) 



/mS^.(iu; + 0+,T) 



2„,2 



327r \ 
AVT 



4M2 (T) X^^T , 
1 T-^^ + — In 



In 



1 _ f-P^^p 



-I3ujp 



UJ 



e(cu2 -4M^(T)) 



(3.40) 





1 4M|(T) 


2 2\ 


^2 _ p2 



(3.41) 



The processes that contribute to the first term above are the creation of two mesons in the 
medium with Boltzmann weight (1 + n^){l + n^^^) minus the inverse process with weight 
ngn^^^, whereas the processes that contribute to the second term are scattering off quanta 
in the bath with weight 



'1 + n^)n^j^^ minus its inverse with weight 



26,16 



Notice that whereas the first term displays the usual two particle threshold, the second 
term is reminiscent of Landau damping in that it only has support below the light cone, for 
\p\ > \uj\ 



10,16 



Neither term contributes on-shell. However there is a new contribution of 
0{X) to mass renormalization and a finite contribution of the same order to wave function 
renormalization. In terms of the renormalized mass, the position of the resonance obtains 
a correction of 0{X) but to this order the width remains the same since the imaginary part 
on shell is already of (9(A^). 



12 



The results of this section also apply if there is a trilinear coupling of the scalar field $ 
to lighter scalars of the form g^x^- this case the result ( p.40| ) applies (with A^f ^ Ag"^) 



but with the mass Mfi{T) in the argument of the square roots in equations ( |3.4CI| , p.41| ) 
replaced by m^(T) the (thermal) mass of the x particle in the loop. If Mr{T) > 2m^(T), 
decay is kinematically allowed (on-shell) and the scalar particle will acquire a contribution 
from the imaginary part of the self-energy on shell from the decay process $ — > XX minus 
the inverse (recombination) process xX ^ ^- The decay (minus recombination) width for 
weak coupling is therefore (see also 



1 - 



4m2(T) g^T 



+ 



In 



-I3W+ 



4m2(T) 



2 2^ Ml{T) 



(3.42) 
(3.43) 



Now the coUisional and the decay width add up. The Landau damping term does not 
contribute to the width of the particle directly, because it has support only below the light 
cone. Notice, however, that the Landau damping contribution to ( |3.42| ) is the same as 
the finite temperature contribution above the two particle threshold in the first term of 
( p.40| ). The sum rule ( |3.27| ) ensures that the Landau damping contribution is "borrowed" 
from the spectral density above the two particle threshold and thus contributes indirectly 
to the width of the particle and to the wave function renormalization (the weight of the 
quasiparticle pole). 

The main purpose of studying the broken symmetry phase is to compare to the kinetic 
description and recognize the terms corresponding to decay and recombination processes in 
the evolution equation for the distribution function to be studied in the next section. 



IV. KINETIC THEORY 

In this section we obtain the evolution equations for the distribution functions of quasi- 
particles, including off-shell effects. Here we present an alternative to the usual derivation 
in which correlation functions are written in terms of relative and "center of mass" space- 
time coordinates and the Wigner transform in the relative coordinates is performed. In this 
approach the Boltzmann equation is obtained in a gradient expansion assuming that the de- 
pendence on "center of mass" coordinates is weak. We will not assume such a situation, but 
instead analyze which are the relevant time scales over which a coarse graining procedure 
must be implemented. 

Let us begin by writing the Hamiltonian of the theory as 

H = Ho + H,nt (4.1) 
Ho = ^J d^x {n^ + (V<l>)2 + (T)$2} , (4.2) 

Where the mass counterterm has been absorbed in the interaction. The interaction part 
Hint will depend on whether T » or T < Tc is considered. The part of the Hamilto- 
nian Hq describes free quasiparticles of renormalized finite temperature mass Mji{T) and 
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is diagonalized in terms of creation and annihilation of free quasiparticle operators ; 
respectively. The pole mass of the quasiparticles will acquire corrections in perturbation 
theory, but these will remain perturbatively small at large temperatures. 

With this definition, the lifetime and weight of the quasiparticles will be a consequence 
of interactions. In this manner, the hard thermal loops which in this theory amount to local 
terms, have all been absorbed in the definition of the thermal mass, which guarantees that 
the microscopic time scales are explicit in the quasiparticle hamiltonian. Such an approach 
has been previously advocated by Lawrie Jl^ . 



Now consider the case in which at an initial time (to be taken ast = to) the density matrix 
is diagonal in the occupation number of the free quasiparticles, but out of equilibrium, with 
initial occupation numbers Nk{to). A comment is in order here. In an initial situation out 
of equilibrium the temperature is not a meaningful quantity and conceptually Mji{T) in the 
above Hamiltonian Hq is not well defined. However we are interested in small departures 
from equilibrium as is implicit in a kinetic description, thus the value of the temperature 
in the resummed mass term should be thought of as the final equilibrium temperature. In 
particular for the relaxation time approximation a typical situation considered is that only 
few modes are out of equilibrium whereas all other are in equilibrium and serve as the "bath" . 
In this case the temperature T in the mass parameter is that of the modes in equilibrium. 
This will become clear in the relaxation time approximation discussed later. 

The Heisenberg field operators at time t are now written as 

$(f , t) = ^ E ^k{t) e*^^ , Mt) = (akit) + al.it) 



where VL is the spatial volume and the time evolution of the creation and annihilation oper- 
ators is through the time evolution operator, uj^ are the same as in ( ^.13| ). 

The expectation value of particle number operator A^^ can be expressed in terms of the 
fields $fc and that of the conjugate momentum H^ as follows, 

N,{t) = (4(t) a,{t)) = -^(Hfc(t) H_,(t) + ul $,(t) $_,(t) ) - \ (4.4) 

where the bracket (• • ■) means an average over the gaussian density matrix defined by the 
initial occupation numbers A^fc(to)- The time-dependent distribution ( [4.4| ) is interpreted as 
the quasiparticle distribution function. 

As in our study of relaxation in the previous section, we examine separately the cases 
T » Tc and T << since the interaction vertices are different (we could study the 
situation in the most general case, but prefer to study in detail the cases separately to 
illustrate clearly the different processes). 



A. T » T, 



c 



The interaction Hamiltonian in this case is given by 
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kl,k2,k3M 



(4.5) 



Taking the derivative of Nk{t) with respect to time and using the Heisenberg field equations, 
we find 



Nk{t) 



1 



2u;k 



^(($)|(t)n_,(t) + n,(t) ($)i,(t)) 



Sm'{T) (<i>fc(t) n_fc(t) + n,(t) $^fe(t)) 



where we use the compact notation: 



mm 



n 



J2 h,+k,+k,-kfl $fci(t) ^kM ^^k.it) 



(4.6) 



(4.7) 



kiMM 



In a perturbative expansion care is needed to handle the canonical momentum (11 = $) 
and the scalar field at the same time because of Schwinger terms. This ambiguity is avoided 
by defining 



{u,{t) mi.it)) = tr p(to) n,(t) (<|.)i,(t) 



lim — — tr 

t^t' d t' 



i^^r-kit)p{to)<^kit')] (4.8) 



Where we used the cyclic property of the trace and the ± superscripts for the fields refer 
to field insertions obtained as variational derivatives with respect to sources in the forward 
time branch (+) and backward time branch (— ) in the non-equilibrium generating functional 

We now use the canonical commutation relation between 11 and $ and define the mass 
counterterm 5rn?{T) = (A/6)A(T) to write the above expression as 



Nk{t) 



X 



.A 



+ 3i 



((*"'Wf.<^>:.(0) + A(T)($+(t)$:,(tO) 

A(T) 



t=t' 



(4.9) 



The right hand side of eq. ( [4.9| ) can be obtained in weak coupling expansion in A. 
Such a perturbative expansion is in terms of the non-equilibrium Green's functions ( |2.4| - 
p.9|) with the basic Green's functions given by (|2.11| , p.l2|) but with the non-equilibrium 
occupation number Nk{to) replacing the equilibrium one n^. An important point to notice 
is that these Green's functions include the proper microscopic scales as the contribution 
of the hard thermal loops have been incorporated by summing the tadpole diagrams. The 
propagators entering in the calculations are the resummed propagators. The terms with 
A(T) are required to cancel the tadpoles to all orders. 

As will be discussed in detail below, such an expansion will be meaningful for times 
t <^ Tj. = \{Nk{t) / Nk{t))\, where is the relaxation time scale for the non-equilibrium 
distribution function. For small enough coupling we expect that will be large enough 
so that there is a wide separation between the microscopic and the relaxation time scales 
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that will warrant such an approximation (see discussion below). At order 0{X) the right 
hand side of (|4.9|) vanishes identically. This is a consequence of the fact that the initial 
density matrix is diagonal in the occupation number. Out of equilibrium, when the equal 
time correlation functions in (|4.9| ) are time dependent, the mass counterterm will acquire a 
(weak) time dependence through the non-equilibrium distribution functions, thus A(T, t) will 
be a slowly varying function of time on a scale determined by r^. Thus from the formidable 
expression ( |4.9| ) only the first term remains after the hard thermal loop resummation and 
we find one of our main results: 



A d 

Qujk df 



t=t' 



(4.10) 



with the understanding that no tadpole diagrams contribute to the above equations as they 
are automatically cancelled by the terms containing A(T) in ( [4.9|) . 

To two-loop order, the diagrams that contribute to (|4.1CI|) are shown in figure 1. We find 
the following expression for the evolution of the distribution function: 



\2 1 

iVfc(t) = — ^ 
^ ' 3 2ujk 



d^ki d^k2 d^k^ 



(27r)3 (27r)3 {2Txf 8uJk,uJk2UJk3 



{27TYS%k -ki-k2- fcg) 



sin [{iJk + ^fci + uJk2 + (^k3)it - to)] 

+ ^fcl + ^^k2 + ^k3 
3 sin [{Uk + UJkl + ^k2 - ^^k3){t 



+ 



^k + ^kl + ^k2 — ^k3 

3 sin [{uJk + ^ki - i^k2 - i^k3){t 



+ 



^k + ^kl — ^k2 — ^k3 

sin [{uJk - iOki - uJk2 - uJk3){t 



+ 



^k ~ ^kl — ^k2 




(4.11) 



where 

M(t) - 

M2{t) - 

Mt) - 
M(t) - 



;i + Nk{t)) (1 + Nki{t)) (1 + Nk2{t)) (1 + Nk3{t)) - Nk{t) Nki{t) Nk2{t) Nk3{t) 
;i + Nk{t)) (1 + Nki{t)) (1 + Nk2{t)) Nksit) - Nkit) Nkiit) Nk2{t) (1 + Nk3{t)) 
;i + Nk{t)) (1 + Nki{t)) Nk2{t) Nk3{t) - Nkit) Nkiit) (1 + Nk2it)) (1 + Nk3it)) 
;i + Nkit)) Nkiit) Nk2it) Nksit) - Nkit) (1 + Nkiit)) (1 + Nk2it)) (1 + Nksit)) (4.12) 



There are several noteworthy features of this expression: 

i) Nkit) vanishes at t = to- This is a consequence of the choice of the initial density matrix 
which is diagonal in the Nkito) basis, so that the occupation number operator commutes 
with the density matrix at t = to- 

ii:) consider times much larger than the time scales implicit in the arguments of the sine 
functions, typically t — to ^ ^R^iT)- Assuming that the rate of change of the occupation 
numbers is very slow so that during this time the occupation numbers have not changed 
significantly, namely, ~ A^fc(to)/iV(to) ^ (t — to) ^ M^^(T), we can approximate 



sin(Vr(t - to)) 
W 



TrSiW) 



(4.13) 
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where W is any of the combinations entering in the arguments of the sine functions above. 
This approximation is the same as that invoked in time-dependent perturbation theory 
leading to Fermi's Golden Rule. Under the assumptions leading to the approximation ( [4.13| ) 
we recognize that each term in the expression (|4.11|) can be identified with a term having the 
same structure in the imaginary part of the self energy given by eq. ( p.21[ ), including off-shell 
processes when the non-equilibrium occupation numbers Nk{to) are replaced by Uk- The first 
term describes the creation of four particles minus the destruction of four particles in the 
plasma, the second and fourth terms describe the creation of three particles and destruction 
of one minus destruction of three and creation of one, the third term is the scattering of 
two particles off two particles and is the usual Boltzmann term. If the Nk{tQ) coincide with 
the equilibrium distribution functions Uk then all the terms in eq. ( [4.11|) vanish identically 
and Nk = 0, including the off-shell terms as a consequence of the energy conservation delta 
functions. 

iii) Using the approximation ( 4.13 ) and assuming that all the modes but the one with 
wavector k are in equilibrium, and the mode with wavevector k is slightly off-equilibrium, 
that is Nki = Uki ; Nk = Uk + Snk{to) to linear order in 6nk{to), we obtain to linear order in 
the Snk{to) the "relaxation time approximation" 

— — = OUkito) (4.14) 

at Uk 

where Ilj{uJk) is given by eq. ( p. 21 ). Thus the relation between the rate of change of the non- 
equilibrium distribution function in the relaxation time approximation and the "damping 
rate" of quasiparticles becomes explicit. 

iv:) Only the third (scattering) term in the evolution equation ( ^.11| ) conserves the total 
number of particles N = Y,k ^k- The other three terms are a consequence of the fact that the 
interaction does not conserve the particle number and this gives rise to off-shell processes that 
change the particle number. However as explained above, these processes only contribute on 
time scales comparable to the microscopic scale. These particle-number-changing processes 
produce departures from chemical equilibrium on short time scales. 

Obviously, equations ( [4.11| , [4.14|) describe only the early time evolution of the distribution 
functions. They neglect the change in the initial occupation numbers and, as they stand, 
cannot be extended to long times. Consider however the following procedure: assume that 
there is a wide separation of time scales in the sense that Tr ^ M^^{T), and consider 
integrating ( [4.11| ) (or the linearized version ( [4.14|) ) in a time interval from t = t^ to tQ + At 
such that M^^(T) ^ At ^ r,.. In this time interval we can approximate the sine functions 
by energy conserving delta functions using ( [4.13| ) keeping the occupation numbers at their 
initial value at t = to- At the end of this interval "reset" all the occupation numbers to the 
values Nk(to + At), and furthermore reduce the density matrix at this time (to + At) to the 
diagonal elements in the basis of the occupation numbers at this time. Then the Green's 
functions will look exactly the same as those given at the initial time t = to but in terms 
of the occupation numbers at t = to + ^t- In this "coarse grained" density matrix there 
are no off-diagonal matrix elements of the creation and destruction quasiparticle operators 
at t = to + ^t- In particular, this reduction of the density matrix neglects new correlations 
of the form (afc(t)a_jt(t)) ; (afc(t)aLfc(t)) which will be generated. Now we can iterate the 
procedure obtaining the equivalent of eq. ([4.11| ) but with the Nk{to + At) on the right hand 
side, iterating this procedure for all times is equivalent to the non-linear equation 
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1 f d^ki d^k2 d^ks 1 ..^N3e3. 



m) = T ^ / 7^ 7^ 7^ ^ i2ny6\k -k,-k,- k,) 

X [6{uJk + UJkl + UJk2 + UJksWlit) 

+3 5{uJk + UJkl + UJk2 - UJk3W2{t) 

+3 6{uJk + UJkl - UJk2 - uJksWsit) 

+S {uJk - UJkl - UJk2 - UJk3)M4{t)] (4. 15) 

Linearizing the above equation, considering only the mode with wavevector k to be 
shghtly out of equihbrium whereas all the other modes are in equilibrium leads to the 
"relaxation time approximation" 

f!^ = -^lii^SnM (4.16) 

dt UJk 



The only term that gives a non-zero contribution in ( [4. 151 ) is the third ("scattering") 



term, because for all the other terms the S functions have no support on the mass-shell. 
Therefore we obtain the Boltzmann equation (|4.15|) and its "relaxation time approxima- 
tion" ( |4.16| ) (where only the second term in ( |3.21| ) contributes to T,i{uJk))- This result 



and the interpretation of the damping rate Tjj{uk)/2ujk as (half) the relaxation rate of the 
quasiparticle distribution function was proposed by Weldon To our knowledge it has 
not been proven before by obtaining the Boltzmann equation directly from a first principle 
calculation and analyzing the "coarse graining" approximations that lead to this result. 

The linearized approximation gives the time scales for relaxation for situations close to 
equilibrium. In the case of small momentum or T >> M^{T) » p we obtain 



MR{T)Tr oc ^ > 1 (4.18) 

where we have assumed very weak coupling. 

The collisional relaxation rate obtained from the hard thermal loop resummation 

^3/2rp 

is very different from the "naive" scattering rate that one would write down for a "classical" 
Boltzmann equation. Such a collision rate would be obtained as 

Tela = (nav) (4.20) 

where n is the number density, cr the 2 — > 2 scattering cross section, v the mean velocity of 
the colliding particles and (■ ■ ■) the ensemble average. For T ^ m one would obtain n ^ T^, 
f fti 1, the Born cross section is (j{E) ^ \^ / E"^ where E is the total energy of the colliding 
particles, which in the high temperature plasma is ~ T. Such a "classical" estimate would 
lead to a scattering rate T^a ~ A^T which for very weak coupling grossly underestimates 



the correct resummed result (|4.19|) for long wavelength modes. 
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B. T <Tc 



In this section we obtain the Boltzmann equation in the broken symmetry case. After a 
shift of the fields on both branches by the vacuum expectation value as v + the 

interaction Hamiltonian contains new terms as compared to the unbroken symmetry case of 
the previous section. These are given by 



1 



fcl,fc2,fc3 



<T)^ $^ ^ L 7ia/Q$- 



h 



6ml 



(4.21) 
(4.22) 



There is now a one loop contribution to the Boltzmann equation. We will not fix the v.e.v. 
to the tree level value. Rather, h is considered as a source term that will be used to cancel 
all the tadpoles arising from the cubic interaction, whereas the mass counterterm will cancel 
the tadpoles corresponding to mass corrections. In this manner, using h to cancel the cubic 
tadpoles the result for v is the true v.e.v. including quantum (and thermal) corrections. 

We now proceed in the same manner as in the unbroken symmetry case of the previous 
section and obtain the evolution equation for the occupation number by using the Heisenberg 
field equations and the non-equilibrium Feynman rules. Using the canonical commutation 
relations between the field and its canonical momentum, and requiring that h cancels all the 
tadpoles arising from the cubic interaction, and 5rn?{T) cancel the tadpoles corresponding 
to mass insertions, we find the new contribution to the evolution equation 



Xv d 



2LJk dt 



7 [((<f+W)^<f 



t=t' 



kl+k2 ,k 



(4.23) 
(4.24) 



fcl,fc2 



This contribution must be added to the obtained in the previous section and given by 
equation ( [4.10|) . To one loop order, the contributions are shown in figure 2 and to this order 
we obtain the following evolution equation for the distribution functions: 



2uJk 



X • 



(27r)3 (27r)3 4cu,icu,2 
sin [{uk + uJki + ujk2)it 



{2'nf5\k-ki 



ko 



+ 



+ ^kl + ^k2 
sin [{Uk - ^fcl - ^k2){t - ^o)] 
^k — ^kl — ^k2 

sin [{ujk + tUfci - ^k2){t - 



+ 



M^{to 

to)] 



^k + ^kl — ^k2 

sin [{uk - ujki + uJk2){t 



+ 



to)]' 



^k — ^kl + ^k2 



Ml{to) 



(4.25) 



where 
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Af^'it) = (1 + Nkit)) (1 + Nkiit)) (1 + Nk2it)) - Nkit) Nkiit) Nk2{t) 

Af^it) = (1 + N,it)) Nk.it) Nk^it) - Nkit) (1 + iV,i(t)) (1 + Nk2it)) 

X^it) = (1 + iV,(t)) (1 + Nkiit)) iVfc2(t) - iV,(t) N,,{t) (1 + iVfc2(t)) 

AC'(t) = (1 + Nkit))Nkiit) (1 + iV,,2(t)) - iVfc(t) (1 + iVM(t)) Nk2it) (4.26) 

Again under the assumption of a separation of time scales and approximating the sine 
functions by energy conserving delta functions using ([4. 131) as before, and considering that all 



the modes but the mode with wave- vector k are in equilibrium with distributions A^A:i(^o) = 
Hki and Nf^lto) = rik + 5nfc(to) we are led to the relaxation time approximation 

^ . -?SMsn,M (4.27) 
cLt UJ}~ 

with Sj(co'fc) given by the imaginary part of the one loop contribution to the self energy, 
equations ( |3.39|) and (|3.40|) . In the broken symmetry case with only one scalar field, the 



imaginary part Yi"j{ujk) vanishes on-shell because of kinematics. However, it is non-zero in 
the case in which the particle $ can decay into lighter scalars x discussed in section III-B 
and below. 

The discussion of the previous case also applies here, in that the above evolution equations 
neglect the change of occupations numbers on the right hand side of the equations, and are 
thus valid only for early times. However, as discussed in the previous case, if there is a 
separation of time scales such that Tj.Mji{T) >> 1 we can integrate the evolution equations 
for times M]^^{T) « [t — to) « r^, approximate the sine functions by energy conserving 
delta functions, and keep the occupation numbers fixed at their values at to- Furthermore 
reducing the density matrix to the diagonal elements in the basis of occupation numbers 
at time t, we can iterate this equation leading to the usual Boltzmann equation. The new 
terms describe the evolution of the distribution function arising from processes of decay and 
recombination in the plasma. Such an iterative procedure leads to the following one-loop 
contribution to the Boltzmann equation 

\^V^ f d^h d^k2 1 ,3^3, 



2ujk J (27r)^ (27r)^ 4c^fcit^fc2 

X6{ujk + CJfei + UJk2)N'l{t) 

+5{uJk - UJkl - UJk2W2i^) 

+6{uJk + UJkl - UJk2W^{^) 

+5{uJk - oJki + u}k2)K{t) . (4.28) 

This one-loop result also applies if the scalar field $ couples trilinearly to lighter scalars Xi 
namely with a vertex g^x^ and we assume that kinematical decay is possible, i.e. Mfj{T) > 
2m^(T). Under the assumption of separation of time scales, which is the case for small 
coupling, justifying the replacement of the sine by delta functions we obtain a Boltzmann 
equation for evolution of the distribution function for the $ field which has the same form 
as ([4.28|) , but with A^f^/2 replaced by 2g'^ and the functions N'^{t) replaced with 



= (1 + Nkit)) (1 + N^.it)) (1 + N^,it)) - Nkit) N^,{t) NUt) 
mt) = (1 + Nkit)) N^.it) NUt) - Nkit) (1 + NUt)) (1 + NUt)) 
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Arm = (1 + N,it)) (1 + N^.it)) NUt) - N,{t) Nlit) (1 + N^,it)) 

Mm = (1 + Nu{t))N^At) (1 + NUt)) - N,(t) (1 + NUt)) NUt) . (4.29) 

Here Nk{t) and N^{t) are the non-equilibrium distribution functions for the quanta of the 



fields $ and x respectively and u!ki,2 = y ^1,2 + "^x^-^)' '^^^ term proportional to J\f2(t) 
describes the decay process of $ — XX minus the "recombination" XX — ^ in the medium. 
The relaxation time approximation for the evolution of Nk{t) = + Snk{t) while the A^^ 
are the equilibrium functions is given by 

^ = -2T^^^^6nk{t) (4.30) 

with r$„^^ the decay minus recombination rate given by equation (p.42| ). 

Thus, after this analysis the nature and validity of the approximations and resummations 
implied in the usual Boltzmann equation (and its linearized version) become clear: 

i) For high temperatures a hard thermal loop resummation to incorporate the proper 
microscopic time scales. Such a resummation is needed to incorporate the relevant time for 
"coarse graining" . 

ii) A wide separation between the microscopic time scale and the relaxation time which, 
as shown by the estimate ( |4.18|) , is warranted for very weak coupling. Such a case cor- 



responds to the coUisional and decay lifetime of the particle being much larger than the 
microscopic time scale given by the inverse mass of the particle in the medium, namely a 
"sharp resonance" condition. This separation of time scales permits the identification of 
intermediate times M^^(T) <^ t <^ t,. such that the "Fermi Golden Rule" approximation 
( [4.13| ) is justified and energy conservation is enforced. 

iii) A reduction of the density matrix to the diagonal elements in terms of the time 
dependent occupation number. This reduction neglects higher order correlations of the 
quasiparticle operators such as {ak{t)ak(t)) and {a\(t)a\(t)) . Therefore the "coarse graining" 
procedure is not only at the level of spatial and temporal averages but also of reducing the 
density matrix, thus reducing the hierarchy of equations of motion. 

Although all of these approximations are implicit in the usual procedure through the 
Wigner transform and the gradient expansion, our method illuminates all the necessary 
steps and the physical justifications at once. 



C. Improvements: 

This method allows us to identify several improvements to the description the kinetics 
of thermalization and relaxation: 

i) The first improvement is to incorporate wave-function renormalization in the definition 
of the quasiparticles in the Hamiltonian as well as coupling constant renormalization. For 
this we write the renormalized fields, canonical momentum and coupling constant as 

$^ = Z^'/'$; n^ = 4/'n; \r = ZIZ^^\ (4.31) 

Notice that 11^ and are canonical conjugate variables. The "free quasiparticle" Hamil- 
tonian is now written in terms of H^, and M^(T) and the difference between the bare 
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and renormalized quantities is accounted for in the interaction part of the Hamiltonian. The 
quasiparticle distribution Nk{t) is defined as in eq. ([4 .41 ) but in terms of the renormahzed 
fields and canonical momenta. Its evolution equation is obtained from the Heisenberg equa- 
tions of motion with the new terms in the interaction Hamiltonian as described previously. 

There will now be new terms in the equations for the evolution of the distribution function 
that will cancel the (divergent) contributions arising from the wave function and coupling 
constant renormalization. These terms will appear at three loops (for coupling constant) and 
beyond but in terms of the quasiparticle propagators. This is precisely the BPHZ program 
of renormalization but here implemented in a non-equilibrium formulation. Clearly such 
an improvement comes at the cost of a complicated structure but this method allows us to 
incorporate such corrections in a systematic manner. 

ii) Rather than assuming a wide separation of time scales and using the approximation 
(|]T|) leading to energy conserving delta functions, one can integrate the evolution equations 
( [4.11| , [4.25| ) numerically in increments t = + At, "reseting" the values of the occupation 
numbers for the next increment to those obtained in the previous iteration. Such an approx- 
imation can be implemented numerically just as in the usual case. However the equation 
for the evolution of the distribution functions will now have off-shell contributions that will 
modify the evolution on time scales comparable to M^^(T). These contributions include 
information on the preparation of the non-equilibrium state and make the derivative of the 
occupation number vanish at t = to as it must (the same feature appears in the derivation 
of Fermi's Golden Rule). 

In the two-loop collisional contribution, the energy- conserving delta functions introduce 
one constraint. Thus a five dimensional k- integral remains, whereas in the "improved" form, 
without energy conservation one must face the full six- dimensional k-integral. Furthermore, 
this "improved" set of evolution equations include terms that change chemical equilibrium 
(do not conserve particle number) on microscopic time scales. These result from the fact 
that the interactions considered in this study do not conserve particle number and may lead 
to new effects in the quasiparticle distribution at long times. Such a possibility should be 
studied further. 

V. CONCLUSIONS AND IMPLICATIONS 

In this article we have presented a new approach, based on non-equilibrium quantum 
field theory, to study linear relaxation of quasiparticles and the kinetic description of the 
evolution of the quasiparticle distribution functions in a scalar field theory. Our method 
incorporates the relevant microscopic time scales through the resummation of hard thermal 
loops in the description of both relaxation and kinetics. The method leads to a system- 
atic and consistent scheme for obtaining both the linearized equations of motion and the 
equations for the distribution functions. It clearly displays the different "coarse graining" 
approximations usually involved in a kinetic description and their range of validity. It allows 
for a systematic and consistent improvement in the kinetic description, including off-shell 
and renormalization effects, as well as processes that change chemical equilibrium on short 
time scales, and is numerically implementable. These improvements in the kinetic equations 
may prove physically relevant if the time scale for relaxational and thermalization processes 
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is not widely different from the microscopic time scales and off-shell processes must be taken 
into account for a more accurate description. 

As a result of our study of both relaxation and kinetics we established a direct proof of the 
relationship first proposed by Weldon between the damping rate of the quasiparticles 
and the relaxation rate of the quasiparticle distribution function in the relaxation time 
approximation after the coarse graining of the kinetic equations. 

We have studied relaxation and kinetics both in the unbroken symmetry phase for T » 
Tc as well as in the ordered, broken symmetry phase for T < Tg. The hard thermal loop 
resummation is very important in the unbroken phase to obtain the proper microscopic time 
scales and to establish a wide separation of time scales that justifies a kinetic description 
in the weak coupling regime. The broken symmetry case provides a setting to incorporate 
processes in which the scalar particle can decay, via a trilinear coupling, to other lighter 
scalars in the theory. We studied the effect of such processes both on relaxation and kinetics, 
obtaining the kinetic equations that include decay and recombination in the medium. We 
also analyzed the effect of processes that only contribute below the light cone that are similar 
to Landau damping. We pointed out that although these processes do not contribute directly 
to relaxation and damping rates, they contribute indirectly through the sum rules obeyed 
by the spectral function. 

The methods developed here are suitable for treating the evolution of the distribution 
functions for the particles produced during the parametric amplification stage in infiationary 
cosmological scenarios plf-p^ and in the supercooling stage of the quark-gluon transition 
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of the expectation value of the scalar fields (mean fields) into the kinetic equations. If there 
is a separation of time scales between the stage of particle production and the onset of ther- 
malization, these equations will lead to a description of particle production via parametric 
amplification and the coUisional processes leading to thermalization and completion of the 
reheating stage. We are currently applying these techniques to that problem, and we expect 
that these methods will be extended to gauge and fermion theories in the near future. 

Some questions that this study poses that are important are the following: is there 
a kinetic description of "critical slowing down" near a phase transition? can a kinetic 
description be extended to temperatures near criticality when long range correlations begin 
to form and relaxational time scales should diverge? 

These questions are relevant for a deeper understanding of the dynamics of phase 
transitions and the non-equilibrium description of thermalization and relaxation of long- 
wavelength fluctuations. We expect to address these and other related questions in the near 
future. 
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Figure Captions: 

Fig. 1 Top: two-loop diagrams contributing to the equation of motion, the dashed 
hne represents an insertion of the non-equihbrium expectation value. Bottom: two-loop 
diagrams contributing to the Boltzmann equation. 

Fig. 2 Top: one-loop diagrams contributing to the equation of motion in the broken 
symmetry phase. Tadpoles have been absorbed in a mass and v.e.v. renormalization. The 
dashed line represents an insertion of the non-equilibrium expectation value. Bottom: one 
loop diagrams contributing to the Boltzmann equation in the broken symmetry phase. 
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